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The confluence of state-of-the-art electronic-structure computations and modern synthetic materials 
growth techniques is proving indispensable in the search for and discovery of new functionalities 
in oxide thin films and heterostructures. Here, we review the recent contributions of electronic- 
structure calculations to predicting, understanding, and discovering new materials physics in thin-film 
perovskite oxides. We show that such calculations can accurately predict both structure and properties 
in advance of film synthesis, thereby guiding the search for materials combinations with specific 
targeted functionalities. In addition, because they can isolate and decouple the effects of various 
parameters which unavoidably occur simultaneously in an experiment - such as epitaxial strain, 
interfacial chemistry and defect profiles - they are able to provide new fundamental knowledge about 
the underlying physics. We conclude by outlining the limitations of current computational techniques, 
as well as some important open questions that we hope will motivate further methodological 
developments in the field. 



Transition metal oxides exhibit the desirable combi- 
nation of high electronic polarizability, originating in 
the chemistry of the transition metal-oxygen bonds, and 
strong electron correlations, from the localized and inter- 
acting transition metal d electrons. As a result of this com- 
bination, the energetics of various interactions - such as 
Coulomb repulsion, strain, orbital bandwidths and Hund's 
exchange - tend to be of similar magnitude. While in 
"conventional" materials, such as semiconductors or metals, 
one of these energy scales dominates and determines the 
macroscopic properties, in transition metal oxides they 
compete, leading to strong lattice-electron, electron-spin, 
and spin-orbit couplings (Figure [lj. The resulting ground 
states tend to have multiple low energy competing phases 
and in turn exhibit enhanced susceptibilities to small 
external perturbations.^ Formation of ABO3 perovskite 
oxides in thin-film form affords an additional parameter 
for controlling the delicate balance among the interac- 
tions to produce unique collective phenomena; indeed, 
drastic changes in properties are reported for thin-film 
oxides, such as the appearance of magnetism in otherwise 
non-magnetic materials^ or the activation of improper 
phase transitions.^ In addition, thin films provide an ap- 
propriate architecture for electric-field-tunable electronic, 
magnetic, and structural phase transitions, and ultimately 
are suitable for technological device integration^ 

Despite the experimental progress in achieving high 
quality coherent perovskite oxide thin films^ and 
heterostructuresj^ there are no general rules for pre- 
dicting the electrical, magnetic, or optical responses at 
oxide heterointerfaces given the known properties of the 
bulk constituents.^ This complication is due in part to 
our limited knowledge of the structure of oxide thin films 
- in particular, oxygen positions are highly non-trivial to 
determine using standard diffraction techniques. In addi- 
tion, the closely competing energy scales which lead to the 
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FIG. 1. The energies of the various interactions listed around 
the circle are all of similar magnitude and the resulting com- 
petition between them leads to strong couplings between the 
electron, spin, lattice and orbital orderings illustrated at the 
vertices of the tetrahedron. The couplings in turn give rise to 
the diverse functionalities of transition metal oxides, such as 
ferroelectricity, colossal magnetoresistance, and superconduc- 
tivity. 

desired novel functionalities, in turn cause the properties 
to be strongly dependent on small changes in atomic struc- 
ture and therefore hard to predict. Thus, while oxide thin 
films have the potential to revolutionize the electronics 
industry through, for example, next generation Mottronic 
devicespES or could provide efficient alternatives for our 
growing energy needsp^EH their adoption in practical 
devices has been slow. This is unlikely to change until 
a detailed microscopic understanding of the atomic and 
electronic structures in oxide thin films is developed. 
This review addresses the role of electronic-structure 
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calculations based on density functional theory (DFT^E^l 
in confronting the complex theoretical challenge posed 
by oxide thin films and heterostructures. Although semi- 
classical phenomenological models and well-developed 
theories for correlated electronic states have been used 
to describe some oxide-oxide interfaces EMBUS] there re- 
mains no consensus as to which of these models (if any) 
is most appropriate for a general description of the elec- 
tronic structure of oxide heterointerfaces. A particular 
deficiency of methods that rely on model Hamiltonians 
is that one of the energy scales in Fig. [T] is assumed to 
be dominant, and the delicate interplay between multiple 
competing interactions is difficult to capture. Also the 
atomic structure of the interface layers must be deter- 
mined (or assumed) for input into the calculation. In 
contrast, DFT-based techniques include all of the quan- 
tum mechanical interactions described by these models 
and the atomic structure on equal footing, provided that a 
suitable exchange-correlation potential is available. Con- 
sequently, they are able to directly explore the funda- 
mental physics^23 of oxide heterointerfaces without a 
priori assumptions about which interactions or struc- 
tural distortions dominate the behavior. As an exam- 
ple, DFT calculations have identified the critical role 
of strain-induced tetragonal distortions in coherent thin 
films that led to dramatic enhancements in heteroepitaxial 
ferroelectricity p * 31 ! superconductivity^ or spin-phonon 
coupling^! depending on the material chemistry. It is 
unlikely that the importance of the tetragonal distortion 
would have been identified if it had been required as 
an input to the calculation rather than obtained as an 
output. 

In this Review, we survey the current capabilities of 
state-of-the-art electronic structure approaches, and re- 
view their application to predicting and understanding 
how strain, coherency and interfacial chemistry combine 
or compete to modify the properties of oxides in thin films 
and heterostructures. We conclude by suggesting future 
research directions and open questions that electronic- 
structure calculations could assist in resolving. 



I. BACKGROUND: STRUCTURAL 
DISTORTIONS IN PEROVSKITE OXIDES 

Before beginning our discussion of the structures of 
thin film perovskite oxides, we briefly review the com- 
mon structural distortions that occur in bulk perovskites. 
The detailed structural distortions that are adopted by 
perovskites are highly significant because they have a 
profound influence on the electronic properties. 

The ideal ABO3 perovskite structure is simple cubic, 
with space group Pm3m [Figure [2J a)]. It consists of 
octahedrally coordinated -B-site cations (usually transi- 
tion metals) with three-dimensionally corner-connected 
BOq oxygen octahedra resulting in • • • O— B— O— B— O- • • 
chains with 180° B—O—B bond angles. Larger cations 
occupy the high symmetry positions of the cuboctahe- 
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FIG. 2. Octahedral rotation phase space in perovskite tran- 
sition metal oxides, (a) Rotations of the octahedra can be 
decomposed about orthogonal axes which intersect at the 
transition metal center, (b) Representative octahedral tilt 
patterns described in the text: a°a°c + and a°a°c- correspond 
to in-phase (c + ) and out-of-phase (c~) rotations of the octahe- 
dra about the z-axis, respectively, while a + a + c° and a~a~c° 
produce similar rotations of the octahedra in the xy-pl&ne. 
The relative rotation of adjacent octahedra from one layer to 
the next is clearly seen for the a°a°c~ tilt pattern (far right). 



dral vacancies between the octahedra (the A-sites). Few 
ABO3 oxides (including the prototype mineral perovksite, 
CaTiOa) in fact adopt this ideal structure, however, and 
in practice, most perovskites exhibit various structural 
distortions that lower the symmetry of the system from 
that of the cubic aristotype. 

The most widely occurring distortions are rotations 
or "tilts" of more-or-less rigid oxygen octahedra around 
one or more high symmetry axes. T hese a re conve- 
niently described using Glazer notation 1 34 * 35 ! in which 
the tilt system is written as a^b^c^ where the letters 
specify the rotations about each pseudo-cubic axis (Fig- 
ure^, and the superscripts indicate whether adjacent 
octahedra rotate in-phase/ferrodistortively (+), out-of- 
phase/antiferrodistortively (— ), or not at all (0). Note 
that the decomposition relies on the octahedral units ap- 
proximately maintaining their regularity, while strictly 
keeping their corner connectivity. As we will see later, 
this picture - which is compatible with Pauling's rulea^ 
for ionic compounds - is approximately correct in most 
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cases. If two letters are the same, then the magnitude of 
the octahedral rotations, regardless of whether they are 
in- or out-of-phase, are equal. A common misconception 
is that the same letters in this notation imply that the cor- 
responding lattice parameters of the crystal are identical, 
which is in fact not the case — it indicates that the nearest 
neighbor transition metal distances along that direction 
are equivalent. Since the octahedra are connected in three 
dimensions, a rotation or tilt in one direction restricts the 
allowed tilts and rotations in other directions. In fact only 
23 tilt systems can be obtained, belonging to 15 unique 
space groups we show examples of some common types 
in Figure [2j and how the combination of tilt systems leads 
to symmetry lowering of the cubic Bravais lattice. 

In addition to these octahedral rotations which 
are driven largely by geometric and electrostatic 
considerationsp^Hil electronically- drive n distortions, par- 
ticularly those caused by the first l 42 l 4 3 l and second-order 
Jahn- Teller effects^HUl are important in determining a 
perovskite's structure. First-order Jahn- Teller distortions 
occur when an electronic degeneracy usually associated 
with the cZ-electrons on the -B-site cation can be removed 
by an appropriate structural distortion. This typically 
manifests as an elongation of some B—O bonds and a 
shortening of others. The associated arrangement of the 
elongations - called a cooperative 3 aim-Teller distortion, 
or an orbital ordering - determines the resulting sym- 
metry of the system. Finally, relative displacements of 
cations and ions that result in polar ferroelectric distor- 
tions further lower the crystal symmetry into a polar 
space group.^ Since these distortions can be described to 
second order in perturbation theory they are often refered 
to as second-order Jahn- Teller effectsPS 



II. STRAIN AND INTERFACE ENGINEERING 
IN THIN FILM PEROVSKITES 

One of the primary routes to engineering the properties 
of a perovskite oxide in a thin film is to leverage the 
elastic strain energy imposed by the constraint that a 
coherently grown film an d its substrate have the same 
in-plane lattice parameters P 4 ^ Appropriate choice of the 
mismatch between the lattice constants of the substrate 
and the film, as well as their relative orientations, can 
be used to impose a specific amount of strain on a film 
by the substrate. While it is widely believed that the 
strain acts by imposing a new in-plane lattice constant 
on the film, exactly how that change in lattice constant 
is accommodated is unclear and difficult to determine 
experimentally. We illustrate this point in Figure [3] One 
possibility, shown in (a) and (b), is that the change in 
in-plane lattice parameter is accommodated entirely by 
a change in the in-plane metal-oxygen bond lengths. In 
panels (c) and (d) we show the other limit: the lattice 
mismatch is accommodated by a change in magnitude (or 
type) of the tilt patterns through rigid rotations of the 
oxygen octahedra, and B—O distances remain unchanged. 



Clearly the two responses will have drastically different 
effects on the functionalities of the film. For example, 
changes in the B—O bond length will affect the magni- 
tude and symmetry of the crystal field splitting, whereas 
changes in B—O—B bond angles determine the strength 
and the sign of magnetic superexchange interactionsPSEU 
Note that, in these simple cartoons the positions of the 
A and B cations are identical in the two limiting cases. 
Since, as we mentioned previously, oxygen positions are 
difficult to determine experimentally, it is difficult to dis- 
tinguish between the two strain-accommodation limits 
from experiments that yield only the cation positionsP^H^ 
We will see later that in most practical cases the actual 
response is intermediate between these two limits. 

Another possible film response to the new substrate- 
enforced lattice parameters is that the film changes its 
equilibrium stoichiometry or defect concentration. Oxy- 
gen vacancies are a particularly common point defect in 
perovskite oxides and it is well established that materials 
with larger concentrations of oxygen vacancies have larger 
lattice constants.^ Since imposition of different strains 
requires growth on different substrates, and associated 
changes in growth parameters, it is once again difficult to 
establish experimentally whether changes in defect con- 
centrations are an intrinsic thermodynamic response to 
strain, or arise from extrinsic factors during processing. 

In addition to the change in the in-plane lattice param- 
eter associated with coherent growth on a substrate, the 
details of the interfacial chemistry and structure are also 
likely to influence the properties of the film. Here possible 
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FIG. 3. In coherently strained perovskite films, the BO& 
octahedra can distort through contraction (a) or elongation 
(b) of the equatorial B-0 bond lengths d due to compressive 
or tensile strain, respectively. Simultaneously or alternatively, 
the octahedra can accommodate the substrate-induced change 
of the in-plane lattice parameters by rotation perpendicular 
to the substrate as in (c), and/or about an axis parallel to the 
substrate plane (d). 
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effects include propagation of a ti lt pattern associated with 
the substrate into the filmj 60 * 61 * chemical bonding across 
the interface] 62 1 63 1 and/or interfacial electrostatics!^^ 
In the following sections, we discuss how first-principles 
calculations can identify the microscopic origins for the 
macroscopic material behavior by decoupling the effects 
of bi-axial strain, symmetry and chemical bonding across 
a perovskite oxide substrate/film interface. 

III. SIMULATION OF THIN FILM EFFECTS - 
HOW DO THE CALCULATIONS ACTUALLY 
WORK? 

In this section we describe the practicalities of how den- 
sity functional calculations for oxide thin films are carried 
out, with a particular focus on their unique capability of 
decoupling the various competing effects that can influ- 
ence a film's behavior. We begin with a brief review of 
the density functional formalism and then show how such 
calculations - through proper choice of simulation cells 
and elastic or electric boundary conditions - can disentan- 
gle the role that epitaxial strain and interface chemistry 
have on the macroscopic properties of perovskite oxide 
heterost ruct ures. 



A. Density Functional Theory 

Formalism. Within the density functional framework, 
the ground state properties of a material are obtained 
through calculation of its electron density, p(r) , which 
uniquely defines the energy, E, of the systemPSl 

E[p(r)] = F[p(r)} +Jdr V ext (r)p(r) . (1) 

Here .F[jo(r)] is a universal functional describing the in- 
ternal quantum mechanical interactions of the electrons, 
and JdrV ex t(r)p(r) is the external potential acting on 
the electrons from the nuclei and any external fields. 

The electronic ground state is found from Eqn. [I] by 
finding the energy density that minimizes the total en- 
ergy. This is formally equivalent to solving the many-body 
Schrodinger equation for the fully interacting electronic 
system. However, an exact analytical form for _F[p(r)] is 
not available and in practice the system is approximated 
by non-interacting electrons experiencing an effective po- 
tential that is formulated to capture the important many- 
electron effectsPSThis effective potential contains the non- 
interacting kinetic energy of the effective single-particle 
states, the classical Coulombic interactions between them, 
and a quantum mechanical "exchange-correlation" en- 
ergy term, V XC1 that approximates the remaining quantum 
mechanical electron-electron interactions. 

Choice of the exact form of V xc is particularly 
important in calculations for transition metal oxide 
heterostructures, because the strongly localized transition 
metal 3d and oxygen 2p electrons, combined with 



localization effects introduced by size quantization, lead 
to explicit and strong electron correlations. For example, 
the widely used local spin density approximation (LSDA) , 
which uses a parameterized forrrP3 of the calculated 
exchange-correlation energy of the uniform electron 
gas, is often inappropriate since transition metal oxides 
show large density variations. Better approximations 
are the generalized gradient approximation (GGA)jSSHHl 
which takes into account variations in the electron 
density through gradient terms, or hybrid functionals, 
which combine a small amount of orbital-dependent 
Hartree-Fock exchange with local or gradient-corrected 
density approximations!22H^ Perhaps the most suc- 
cessful simple extension to the LSDA is the LDA+Z7 
method, in which an orbital-dependent energy cost for 
adding an additional electron to an already occupied 
manifold mimics in spirit the Hubbard U Coulombic 
repulsion.^ In all cases, care must be taken that the 
choice of functional is appropriate for the materials to be 
studied if physically meaningful results are to be obtained. 

Extracting material properties. By applying suit- 
able optimization techniques to Eq. [T] the ground state 
charge density and total energy can be readily obtained 
for a fixed crystal structure with a specific lattice geom- 
etry and atomic positions. All ground state properties 
of the system - such as magnetic ordering, densities of 
electronic states, ferroelectric polarizations, etc. - can 
then be directly obtained from the charge density, in some 
cases with external fields includ ed expli citly in the Hamil- 
tonian during the optimization! 6 -^*!*^! j n particular, the 
atomic forces and stresses can be computed, and the ionic 
positions and cell parameters adjusted so that the forces 
and stresses are reduced to zero. This "relaxation" yields 
the lowest energy, ground-state structure. Since in gen- 
eral the relaxation is a complex, multi-variable problem 
with many possible local minima,^ it is usually achieved 
in practice by comparing trial structures with different 
crystal symmetries, obtained by freezing in combinations 
of the unstable phonons calculated for a high symme- 
try reference phase. The lattice phonon frequencies and 
eigenvectors can in turn be used in the interpretation 
of Raman and infrared spectroscopies as well as for the 
study of structural phase transitions. 



B. Homoepitaxial strain. 

We call the first type of simulation that is commonly 
employed the homoepitaxial strain approach. In homoepi- 
taxial strain calculations, we, in fact, model a crystal 
using a simulation cell that is periodic in all three di- 
mensions and subject it to a strain parallel to a chosen 
lattice plane by varying the lattice parameters away from 
their equilibrium values in that plane [Figure [4^ a)]. Such 
simulations then allow the intrinsic role of epitaxial strain 
imposed by lattice matching with a substrate to be deter- 
mined separately from any other effects associated with 
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FIG. 4. Illustration of the homoepitaxial and heteroepitaxial strain models. Possible choices for the fundamental unit cell used 
to represent the homoepitaxial film in a periodic boundary condition DFT-calculation are shown in (b). The left cartoon shows 
a primitive 5-atom perovskite unit cell, and \/2 x V2 x 2 and 2x2x2 supercells are shown center and right. The particular 
cell size is selected in order to accommodate arrangements of various structural and electronic internal degrees freedom. In the 
heteroepitaxial strain calculations, the substrate and film are both present in the primitive heterostructure cell (d). 



the presence of the interface; this is clearly unfeasible 
experimentally. Often bi-axial strain is applied, and in 
this case the strain is given by 



a 

where a is the in-plane lattice parameter imposed on the 
fundamental unit cell [Figure ^b)} in the homoepitaxial 
strain calculation, and ao is the calculated equilibrium 
lattice parameter of the bulk material. Different 
compressive (a < ao) and tensile (a > ao) strain states 
are obtained by varying a. Uniaxial strain can be applied 
by varying only one of the in-plane lattice parameters, 
and anisotropic biaxial strain by varying the two in-plane 
parameters by different amounts. In all cases the 
out-of-plane c lattice parameter and atomic positions are 
relaxed to minimize the stress on the unit cell and forces 
on the ions, respectively. 

Choice of simulation cell. Care should be taken to 
select a simulation cell that allows exploration of all likely 
tilt, rotation, and orbital ordering patterns that might be 
induced by the epitaxial strain constraints. As we will see 
below in the results section, the tilt pattern in the film is 
often different from that in the bulk parent phase, and the 
default primitive simulation cell that correctly accounts 
for the bulk system's structure might not have sufficient 
flexibility to accommodate the lowest energy structure in 
the film. In such cases, a supercell is constructed from 



multiplication (and at times rotation) of the primitive 
cell's lattice vectors. Such examples of supercells that 
can be used to perform homoepitaxial strain calculations 
are shown in Fig. Qb). Each of these simulation cells 
are an integer multiple of formula units (f.u.) of the 
primitive 5-atom perovskite building block. A drawback 
with the larger cells, however, is that they are more 
computationally expensive. 

The importance of selecting the size of the supercell 
is best illustrated through a simple example. The bulk 
low-temperature structure of SrTiOs has an a°a°c~ tilt 
system, which requires a supercell with a \[2 x y2 x 2 
(20 atom) larger periodicity than that of the ideal 5-atom 
perovskite [Fig. Qb, center)]. If epitaxial or uni-axial 
strain were to stabilize the mixed-tilt a°b + c~ pattern 
that is found in SrZr03, however, a 2x2x2 (40 atom) 
supercell would be needed. The smaller unit cell would 
fail to find the global ground state structure and the 
results could be misleading. We survey in the Results 
section, how simulation cell choice is crucial to identifying 
strain-induced changes to the oxygen octahedral rotations 
in perovskites. 

Symmetry consequences. Even at the homoepitax- 
ial strain level, without explicit inclusion of the substrate 
in the calculation, the imposition of specific in-plane lat- 
tice parameters often lowers the lattice symmetry of the 
film. Here, we look briefly at why this occurs and discuss 
the implications of the symmetry change on both the 
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FIG. 5. (a) Structure of a rhombohedral perovskite with 
symmetry R3c exhibiting a~a~a~ octahedral rotations about 
the three-fold [111] axis (red arrow), (b) Projection of the 
three dimensional structure shown in (a) into the ab epitaxial 
plane. In (c) we show the same projection without the oxygen 
atoms for clarity. The blue square shows how forcing the 
rhombohedral cell to match with a cubic substrate reduces the 
symmetry in this case to monoclinic. The effect of different 
strain symmetries in the perovskite film on the 3d one-electron 
orbital level splittings is shown in (d). 



practicalities of the calculations and the macroscopic film 
properties. 

We first describe the possible symmetry modifications 
with the simplest possible example: A material which in 
the bulk has cubic symmetry (space group Pm3m), and 
which is grown on a (001) interfacial plane of a substrate 
also with cubic symmetry. An example of such a material 
would be an ideal cubic perovskite, with no symmetry- 
lowering rotations of the oxygen octahedra, Jahn- Teller 
or ferroelectric distortions. When biaxial strain is im- 
posed (a ^ ao), it is clear that the space group symme- 
try of the material is immediately reduced to tetragonal 
(P4/mmm), with the three-fold axes along the body di- 
agonal of the unit cell removed, but the four- fold and 
mirror symmetry elements preserved. Note that of course, 
the symmetry can be further lowered due to internal 
structural distortions. For example, a common effect of 
in-plane compressive strain on perovskite dielectrics is to 
induce a ferroelectric polarization along the c-axis.^l In 



TABLE I. Symmetry consequences on the film lattice type due 
to two-dimensional thin film epitaxy. In all cases the orienta- 
tion of the substate is given with reference to the pseudocubic 
lattice vectors of the aristotype perovskite phase. Note that 
an appropriately oriented tetragonal substrate will result in 
the same symmetry lowering as a cubic substrate with bi-axial 
strain. 



Substrate 


Parent film under bi-axial strain 


cubic 


rhombohedral 


orthorhombic 


cubic (001) 


tetragonal 


monoclinic 


orthorhombic 


cubic (110) 


orthorhombic 


triclinic 


monoclinic 


rhombohedral (001) 


monoclinic 


triclinic 


triclinic 


orthorhombic (001) 


orthorhombic 


triclinic 


orthorhombic 



this latter case, an additional mirror plane is lost and the 
space group symmetry is reduced further to PAmm. 

We next illustrate a more complex situation in which a 
rhombohedral perovskite - characterized by octahedral 
rotations around the three fold axis along the pseudo- 
cubic [11 Indirection [Figure [5]ja)] - is placed on an (001)- 
oriented cubic substrate. The anti-phase rotation of the 
B0 6 octahedra about the [111] axis leads to a shear dis- 
tortion in the interaxial cell angles. This can be seen 
in Figure [5^b) , in which we have projected the three- 
dimensional structure onto the two-dimensional a6-plane; 
the angle between the a and b lattice vectors clearly 
deviates from 90°. Imposition of a coherent epitaxial 
constraint onto a cubic (100) substrate, however, forces 
the in- plane inter-axial angle to be 90°, as shown in 
Figure [5^c) . Since the epitaxial constraint does not im- 
pose restrictions on the out-of-plane angles, the formerly 
three-fold axis is reduced to a — operation resulting in 
monoclinic symmetry.^ In Table]]] we list the symmetry 
reductions that occur for other common lattice types on a 
range of common substrate symmetries and orientations. 

The reductions in lattice symmetry discussed above of 
course have consequences on the properties. Changes in 
the point symmetry of the transition metal cation site 
can modify the crystal field splitting of the d orbitals. 
For example, in Figure [5^d) we show the effect on 
the d-orbital energy levels of reducing an octahedrally 
coordinated transition metal ion's symmetry to tetragonal 
(left) or trigonal (right). The former is often believed to 
be the experimental situation in ultra-thin perovskite 
films on a square substrate in which pseudomorphic 
growth is maintained.^ In contrast the latter is likely 
for thicker films of perovskites with bulk rhombohedral 
phases in which lattice relaxations have occurred.^ Such 
changes in the crystal-field splitting can have dramatic 
effects in d n perovskite oxides by modfying orbital 
degeneracies and in turn Jahn- Teller and orbital ordering 
patterns, as well as allowing transition metal d-electrons 
with different orbital angular momentum to mix more 
(or less) strongly. 

Finally, we point out that homoepitaxial strain sim- 
ulations offer an exciting complement to experimental 
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studies of strain effects on thin film properties since the 
substrate's lattice parameter a can be tuned continuously 
in the calculations to explore the full strain-structure- 
property phase space. As a result it is possible to identify 
critical strain values at which phase transitions should 
occur and susceptibilities diverge. In contrast, the ex- 
perimental situation tends to be less flexible since the 
choice of substrates on which to grow thin films is limited 
by the availability of high purity single crystals, and the 
equilibrium volumes for their particular chemical composi- 
tions. Of course homoepitaxial calculations omit explicit 
structural and chemical interactions with the substrate 
material, which in some cases might be important and 
even dominate the strain-induced behavior. We discuss 
how such factors are incorporated in so-called heteroepi- 
taxial simulations next. 



C. Heteroepitaxial strain simulations 

In heteroepitaxial strain simulations a two-component 
supercell is used, with a second material included explic- 
itly to model the presence of a substrate [Figure [4^c) and 
(d)]. With both constituents included, it is possible to 
identify how the physical and electronic structure of the 
atoms in the different layers of the film (i), interface (ii), 
and/or substrate (Hi) contribute to the properties of the 
heterostructure. Then, by comparing the results of ho- 
moepitaxial strain calculations with these heteroepitaxial 
strain calculations, the role of the substrate and interface 
can be isolated. 

The first step in a heteroepitaxial strain simulation is 
the selection of a unit cell that is appropriate for study- 
ing the properties of interest. An example is shown in 
Figure |4|d) . Since periodic boundary conditions are used, 
in practice the calculations model a periodic superlattice 
with an infinite array of interfaces rather than a single 
heterointerface. If this is indeed the experimental sit- 
uation, then the repeat unit should be chosen to most 
closely match the periodicity of the experimental super- 
lattice within the limits of available computer resources. 
In this case, it is usually appropriate to relax all lattice 
parameters and internal atomic positions to their lowest 
energy values. If instead the calculations aim to answer 
questions about a single component film on a thick sub- 
strate, the supercell should be chosen to be as large as 
possible so that interactions between the interfaces are 
minimized. In this limit, the in-plane lattice constant 
is usually constrained to that of the substrate, and the 
atoms in the middle layers of the substrate material are 
often fixed to their bulk positions. 

A key question to be addressed in heteroepitaxial calcu- 
lations is whether the tilt and rotation patterns found in 
the substrate template across the interface into the film 
material, and if so, how do these distortions modify the 
macroscopic properties. In Table \H\ we list the octahedral 
tilt patterns adopted by many widely used substrate ma- 
terials at various temperatures. The propagation of these 




FIG. 6. Example of a frustrated antiferrodistortive octahedral 
rotation system in an odd period superlattice (a). The sense 
of the octahedral rotation is indicated by the red arrows, and 
due to the chosen periodicity, a rotation direction cannot be 
selected in the blue block that is compatible with the periodic 
boundary conditions. To eliminate the artificial octahedral 
frustration in the superlattice, a supercell doubled along the 
z-direction is needed to be fully compatible with the a°a°c- 
rotation pattern. 



tilt patterns across interfaces has sometimes been invok ed 
in the literature to explain observed phenomen d 53 * 54 ! de- 
spite being difficult to confirm experimentally. Density 
functional calculations with full or constrained structural 
optimizations are a powerful tool for testing the veracity 
of this assumption. In addition, large superlattices make 
it possible to analyze layer-by-layer changes in the atomic 
and electronic structure to identify whether there is a 
critical thicknesses at which the film recovers its bulk tilt 
pattern and other properties. 

An additional point to note in heteroepitaxial strain 
calculations is that the internal distortions of the atoms, 
primarily the rotations of the octahedra, must also 
be compatible with the size of the unit cell used to 
simulate the superlattice. As discussed in the case of 
homoepitaxial strain, the chosen in-plane periodicity 
must allow the necessary tilt and rotation patterns; the 
supercell shown in Figure [4^d), with its single perovskite 
unit cell in plane, clearly prohibits this. In addition the 
out-of-plane periodicity must be chosen to avoid artificial 
frustration of the tilt system (Figure pi. For example, the 
out-of-phase a°a°c~ tilt pattern requires an even number 
of 5-atom perovskite blocks along the rotation axis to 
accommodate the full periodicity of the zone-boundary 
phonon mode. If, for example, 1/2-period perovskite 
superlattice with only three 5-atom blocks is simulated 
in this tilt pattern, as shown schematically in Figure 

a), the octahedral rotation in the last block becomes 
frustrated: The unit cell below it expects it to rotate 
clockwise, whereas that above it expects it to rotate anti- 



TABLE II. Lattice parameters, crystal structures and tilt systems of common substrate materials used in oxide thin film growth. 
Pseudo-cubic lattice parameters are given in parentheses. For rhombohedral space groups a pc w a/y2. For orthorhombic 
substrates values in parentheses correspond to the average pseudocubic spacing (y/a? + b 2 /2) along the [110] direction. 



Substrate 


Structure 


Temperature 


Tilt System 


Lattice constants (A) 




Reference 


SrTiOg (STO) 


cubic (221, Pm3m) 


> 105 K 


a°a°a 


a = 3.905 




1551 




tetragonal (140, M/mcm) 


< 105 K 


a°a°c- 


c/a = 1.0056 






LaA10 3 (LAO) 


cubic (221, Pm3m) 


> 800 K 


a a°a 


a = 3.81 




|8j 




rhombohedral (167, i?3c) 


< 800 K 


a~ a~ a~ 


a = 5.36 (3.79) 




EZl 


LSAT 


cubic (221, PmZm) 


> 150 K 


a a a 


a = 3.87 




EFJ 




tetragonal (140, M/mcm) 


< 150 K 


a°a°c- 


a = 5.46, c= 7.73 






LaGa03 


rhombohedral (167, R3c) 


> 420 K 


a~ a~ a~ 


a = 5.58 (3.94) 




[89] 




orthorhombic (62, Pnma) 


< 420 K 


a+b-b- 


a = 5.49, b= 5.53, c = 


7.78 (3.89) 


mm 


DyScOa (DSO) 


orthorhombic (62, Pnma) 




a+b-b- 


o = 5.44, 6 = 5.71, c = 


7.89 (3.94) 


ED 



clockwise. A possible solution is to double the primitive 
heteroepitaxial cell along the rotation axis and then per- 
form calculations on the supercell depicted in Figure K3[b). 



IV. RESULTS: HOMOEPITAXIAL STRAIN 

A. "Simple" examples: Strain— octahedral tilt 
coupling in rhombohedral LaAlC>3 and LaNiC>3 

We begin our results section with a discussion of the 
behavior of two ostensibly simple examples: LaAlC>3 and 
LaNiC>3 films strained on (OOl)-oriented cubic substrates. 
Both compounds crystallize in the rhombohedral i?3c 
space group with the a~a~a~ tilt pattern, and in both 
cases, the bulk ground state structure consists of only 
rotational distortions from the ideal cubic structure. The 
bulk pseudocubic lattice parameters are 3.79 and 3.84 A 
respectively, resulting in Al-O-Al (Ni-O-Ni) bond angles 
of 171° (165°) and Al-0 (Ni-O) bond lengths of 1.90 (1.94) 
A. In this symmetry class, the B0 6 octahedral rotation 
angle is known to vary strongly with the shape of the 
Bravais latticef^HHIl Therefore, we expect that changes in 
the degree of monoclinicity introduced by the substrate 
strain (Table [I]) will have a large effect on the internal 
atomic positions. Finally, LaAlC>3 is a robust wide band- 
gap insulator with no tendency to ferroelectric distortion, 
whereas LaNiOa is metallic, with the additional complex- 
ity of partially filled d orbitals, suggesting proximity to 
Jahn-Teller or charge disproportionation instabilities, as 
well as possible metal-insulator transitions! 95 " 97 * These 
factors combine to make LaAlC>3 and LaNiC>3 model sys- 
tems for studying coupling between strain and octahedral 
rotations, and the influence of additional complexities on 
this coupling. 

To investigate the possible rotation patterns present 
in a homoepitaxial strained film, it is useful to first use 
symmetry analysis to determine the possible octahedral 
tilt systems and space groups that are compatible with 



symmetry of the cubic (001) substrate. In this case these 
are: aT aT cr (15 C2/c or 14 P2i/c), a~ aT ~c + (62 Pnma), 
a + a + c~ (137 PAijnmc), a~a°c° (69 Fmmm), a~a~c° 
(74 Imma), a Q a°c° (221 Pm3m), a+a Q c° (55 Pbam), 
a + a + c° (139 IA/mmm), a°a a c" (140 IA/mcm), a°a°c + 
(127 PA/mbm), and a + a + c + (71 Immm). If experimen- 
tal pressure or temperature phase diagrams are available, 
one can also further narrow this set to likely candidates. 
The relative energies of these structural phases are then 
calculated and compared to determine the lowest energy 
structure at each strain state. 

Following identification of the important tilt patterns, 
the change in magnitude of the octahedral rotations about 
each direction for a given Glazer tilt system is commonly 
described in terms of octahedral "rotations" and "tilts." 
These octahedral connectivity descriptors are defined rela- 
tive to the substrate orientation (Fig. [7| : The tilt angle is 
given as (180 — 4>)/2 and the rotation angle is (90 — 0)/2. 




FIG. 7. Definitions for the rotation (~ (j>) and tilt (~ 6*) angles 
that are often used in the literature to describe substrate- 
induced changes to the Glazer octahedral tilt patterns. 
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strain (%) strain (%) 

FIG. 8. Zero kelvin strain-octahedral phase diagrams for two bulk rhombohedral perovskites. (a) Change in the A106 octahedral 
rotation angles in LaAlOs and (b) the energy stabilization for each phase with the corresponding tilt patterns, (c) In LaNi03, 
the a~a~c~ tilt pattern is found to be the ground state for all strains investigated. The "sense" of the a~a~c~ octahedral 
rotation pattern in LaNiOa changes abruptly as the axial ratio approaches unity (d). Data reproduced with permission from 
Refs. 1981 and 1571 Copyright 2010, American Physical Society. 



In this way, the effects of substrate-strain on the octahe- 
dral B-O-B bond angles can be linked to the material 
properties. For a film grown on a substrate with a square 
net (bi- axial strain), the geometric relationships between 
Glazer's rotation angle definitions about the pseudo-cubic 
lattice vectors (Figure [2]) to the rotation 9 and tilt <j) 
angles are as follows: 9 — 7 and <j> — (a + (3)/\/2. 

Using the octahedral symmetry guided approach, Hatt 
and Spaldin mapped out the octahedral-strain phase dia- 
gram for LaAlC>3 [Figure |8ja)], which we reproduce and 
describe hereP^We see in Figure [8ja) that for very small 
strains, the rhombohedral-like (a~a~a~) pattern of the 
bulk system persists, albeit with a monoclinic lattice dis- 
tortion from constraining the in-plane lattice parameters 
and angle. Small compressive strains of 0.2%, however, 
are able to fully suppress the bulk octahedral a~a~a~ 
tilt pattern and change it to a°a°c~, which has rotations 
around the out-of-plane axis only. In this region, the 
lattice is tetragonal with 90° interaxial angles and space 
group I 4/ mem. This behavior can be naively understood 
in terms of the reduction in area perpendicular to the ro- 
tation axis in response to the in-plane compressive strain. 
In contrast, under tensile strain, the a~a~c° pattern is 
stabilized, with rotations around the [110] axis only (space 
group Imma). This rotation pattern reduces the length 
of the film's unit cell along the out-of-plane direction (in 
accordance with elastic theory due to the tensile strain 
in plane) while keeping the Al-0 bonds nearly constant. 

In Figure [8](b) we show the calculated total energies of 
LaAlC>3 as a function of in-plane biaxial strain. We see 
that the energy differences between different patterns of 
tilt and rotation at typical strain values are only a few 
meV/f.u. In Ref. HH]it was also found that the energy of 



unstrained but coherent LaAlC>3 is higher than that of 
the relaxed bulk rhombohedral structure. If all interaxial 
angles are fixed to be 90° - as is often assumed in the 
literature - the increase in energy is 0.6 meV/f.u., whereas 
when a relaxation to monoclinic symmetry is allowed it 
is only 0.2 meV per formula unit higher in energy. This 
suggests that the films will in practice exhibit monoclinic 
distortions. 

Since the nature of the octahedral rotations depends 
on a delicate balance between bonding and electrostatic 
interactions in the solid, we might expect the behavior 
in metallic LaNiC>3 to differ substantially from that in 
insulating LaAlC>3. Therefore, we next describe the octa- 
hedral structure evolution with strain for LaNi03 [Figure 
|8](c)] . While the variations in tilt and rotation angles 
show some similarity to the LaAlC>3 case - the amount of 
[001] rotation is increased with compressive strain and de- 
creased with tensile strain, with the opposite behavior for 
the [110] tilt - in this case the monoclinic C2/c structure 
persists over the whole strain range. A recent detailed x- 
ray diffraction study of strained LaNi03 films on SrTiC>3 
and LaAlC>3 substrates by May and co- workers^ yielded 
refined structures in excellent argeement with these com- 
putational results: neither the tilts nor the rotations 
are ever completely de-stabilized by the substrate-induce 
strain. 

It is clear from Fig. [8] that, while there is no change in 
symmetry for LaNi03 over the strain range investigated, 
there is a distinct structural transition characterized by 
a discontinuity in the magnitude of the tilt and rotation 
angles; this corresponds to an abrupt reorientation in 
the axis about which the NiOg octahedra rotate. This 
first-order phase transition is classified as isosymmetric 
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strain (%) 

FIG. 9. The (OOl)-strain dependence of the in-plane (circles) 
and out-of-plane (diamonds) bond angles and lengths for a 
homoepitaxial LaNiC>3 film. 



control of both parameters, therefore, would be highly 
desirable in attempting to engineer specific behaviors. To 
address this question, we show in Figure [9] the change 
in in-plane and out-of-planc Ni— O— Ni bond angles and 
Ni— O bond lengths as a function of strain. The first 
striking result is that the in-plane bond angle is only 
weakly sensitive to epitaxial strain (0.62°/ percent strain, 
with a discontinuity at the isosymmetric phase transition) , 
while the out-of-plane angle can be tuned by 3.6°/ percent 
strain. Conversely, the in-plane Ni— O bond lengths are 
strongly strain dependent, since they are accommodating 
the change in in-plane unit cell area with strain. The out- 
of-plane bond lengths are only weakly strain dependent 
because changes in the out-of-plane angle take up the 
change in out-of-plane lattice parameter. It is clear that a 
simple picture of strain accommodation via rigid octahe- 
dral rotations is not appropriate, and in fact changes in 
bond lengths — in this case specifically the in-plane bond 
lengths — mediate a substantial portion of the change in 
lattice parameters. 



since the atomic structure remains monoclinic and the full 
symmetry and Wyckoff positions of the C2/c space group 
are retainedPS (We note the authors of Ref. |S7] predict an 
additional phase of lower symmetry (P2\/c) in LaNiC>3 
films under tensile strain which shows a small charge 
disproportionation (CDP) but with negligible changes in 
the rotation angles. It is approximately 2 meV lower in 
energy and would result in a "normal" , non-isosymmetric 
transition.) Interestingly, the isosymmetric phase transi- 
tion occurs where the axial ratio of the crystal approaches 
unity [Figure |8jd)]. This can be understood on geometric 
grounds: when the ideal NiOg octahedra are recovered 
near 0.5% strain, the bi-axial lattice distortion imposed 
on the film by the substrate leads to a tiling of corner- 
connected octahedra that is incompatible with unit cell 
size.™ 

Finally, we note that if the calculations for strained 
LaNiC>3 do not allow for full relaxation of the out-of-plane 
lattice parameter, LaNiOs appears to have the same be- 
havior as LaA103. This is important for two reasons: 
First, it indicates that full structural optimizations are 
essential, and qualitatively incorrect behavior can be ob- 
tained by artificially neglecting some structural degrees of 
freedom. And second, it points to the strong sensitivity 
of the octahedral rotations to the details of the elastic 
response of the material. 

These two examples show that strain clearly couples to 
the magnitude of the octahedral rotations about differ- 
ent crystallographic axes relative to the interface. In the 
next section, we look at how strain influences separately 
the B—O—B bond angles and B— O bond lengths. This 
separation is important because the effects of bond angle 
and bond length changes on electronic properties, such 
as band widths and band gaps, and magnetic properties 
(exchange interactions) are often quite different. Separate 



B. Orthorhombic SrRu0 3 and CaTi0 3 

The "simple" rhombohedral perovskites of the previous 
section were useful examples for illustrating the conse- 
quences of symmetry lowering and elastic energy accom- 
modation in thin films. Here we discuss two examples 
- SrRuOs and CaTiC>3 - which have the more complex 
orthorhombic Prima symmetry, with the a~a~c + tilt 
pattern. This class is particularly important since it is 
adopted by the majority of ABO3 perovskites ) 39 l 10 1 1 The 
thin film behavior of materials with this tilt system is more 
complicated since the structure already has orientational 
anisotropy in the bulk, with the in-phase c + rotations 
along the long axis of the unit cell. As a result the films 
can have two unique orientations on a substrateP^I When 
the long axis lies in the epitaxial plane the film is described 
as a&-oriented and when it is out-of-plane as c-oriented. 
While experimental t hin film s often show a mixture of 
these two orientations LloaiHSl] first-principles calculations 
can determine which is energetically more favorable at 
each strain state, as well as the separate properties of the 
two orientations. In addition, we note that the orthorhom- 
bic structure is reached from the cubic aristotype by the 
softening of two zone-boundary instabilities of different 
symmetry: one at the zone corner (i?-point) and one at 
the edge (M-point). In contrast, the bulk R3c structures 
discussed earlier are reached with only one i?-point insta- 
bility. We therefore expect that the a~a~c + tilt pattern 
might show a quite different strain response — possibly 
with less strong coupling between the rotations and the 
strain. 

We illustrate the effect of strain on orthorhombic per- 
ovskites using two test m ater ials: metallic SrRuC>3, a 
common electrode materiaP^ used in thin film growth, 
and the insulating, prototypical perovskite, CaTiC^. As 
in the previous section, one of our choices is a corre- 
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lated me tal,E IZhi09| anr \ the other is a wide-band gap 
insulatorp^U although we will see in this case that the 
dielectrically active Ti 4+ ion in CaTiOa leads to further 
complexity compared with the inert Al 3+ ion in LaA103. 
SrRu03 (CaTiOs) has rotation and tilt angles of 5.8° and 
7.6° (8.7° and 11.7°), respecti vely, wit h metal-oxygen- 
metal bond angles of 163° (i56°).liillii2lNote, these angles 
are further away from 180° than our rhombohedral exam- 
ples, thus indicating that the compounds are more highly 
distorted. 



In Figure 10 we show the results of h omo epit axial 
strain calculations by Zayak and co-worker:P^ for differ- 
ent orientations of orthorhombic SrRuOa. We see that 
the a~a - c + tilt pattern is maintained for all strain val- 
ues explored. For the c-oriented structural variant, the 
in-plane lattice parameters of the orthorhombic film are 
those that set the a~a~ component of the tilt pattern. 
Since they are already equal in the bulk we expect similar 
behavior to the rhombohedral case. Indeed, as strain is 
increased from compressive to tensile, the tilt (rotation) 
angles increase (decrease) as observed previously. For 
the afr-orientated structure, however, where the lattice 
parameters of the c + rotation axis and one of the a~ 
axes are set by the substrate, a markedly different strain 
dependence of the rotation and tilt angle occurs. In this 
case both tilt and rotation angles decrease (increase) with 
compressive (tensile) strain. This is the first example that 
we have seen of strain giving an overall change in the 
magnitude of the angles; in the previous cases a reduction 
in rotation was always compensated for by an increase in 
tilts. 

Controlled growth of (llO)-oriented orthorhombic per- 
ovskites could therefore provide a route for controlling the 
B—O—B bond angles and as a result the electronic band- 
width in functional oxides. Evidence of such orientation- 
dependent electronic properties is reported for isostruc- 



[110] tilt 




[001] 
rotation 

_j I I 

-3 -1.5 1.5 3 
strain (%) 

FIG. 10. Changes in the RuOe octahedra rotation angles 
along the [001]- (circles) and [110]- (squares) directions for a 
homoepitaxially strained SrRuOs film. We show the strain 
dependence for the orthorhombic films with the long axis 
set perpendicular (c-oriented, filled symbols) and parallel (ab- 
oriented, open symbols) to the (OOl)-epitaxial plane. Data 
reproduced with permission from Ref. 11131 Copyright 2006, 
American Physical Society. 



tural LaTiC>3 thin films here, LaTiC>3 films under 
nominally similar magnitudes of strain on (OOl)-oriented 
substrates show robust metallic behavior, while those 
grown on (llO)-oriented surfaces are highly insulating. 

Zayak and co-workers also calculated the magnetostruc- 
tural coupling in SrRuC>3 and found it to be both sub- 
strate orientation (0.35 /J,s/f.u. difference between (110)- 
and (OOl)-oriented) and strain dependent (0.13 /is/f.u. 
per percent strain) Later they predicted an un- 
usual low-spin S = to high-spin S = 1 state tran- 
sition for the Ru + cation under bi-axial elastic strain 
conditions. 115 While such an on/off control of magnetism 
with strain in strontium ruthenate remains to be con- 
firmed experimentallyp^the prediction is consistent with 
isovalent chemical substitution studies (that mimic the 
application of pressure), which show that the ferromag- 
netic ground st ate is highly susceptible to changes in the 
rotation anglesP^H 

Next we describe the calculated strain behavior of or- 
thorhombic Pnma C aTiC ^. A detailed first-principles 
study by Eklund et a/P^S found similar evolution in the 
octahedral tilts and rotations as those shown in Fig.[T0]for 
SrRuC>3. In addition, however, the authors of Ref. Ill 71 
found that the relative stability of ab- and c- orientated 
structures can be tuned with bi-axial strain. Under com- 
pressive and small tensile strains the a6-oriented films are 
more stable than the c-oriented films. However with in- 
creasing tensile strain, the a6-oriented films are relatively 
de-stabilized with c-oriented films becoming energetically 
more favorable at around +1.5% strain. Intriguingly, at 
very large tensile strains (around +4 %), a ferroelectric 
ground state is obtained in the stable c-oriented films with 
the direction of polarization along the [110]-direction. In 
contrast, no polar phases were found in the a6-oriented 
films between -3% compressive strain through the range 
of stability to 1.5 % tensile strain. Because the change 
in the magnitude of the octahedral rotations under such 
large tensile strains is calculated to be small with respect 
to the strain-free case, the authors attributed the activa- 
tion of the polar instability to a large strain-polarization 
coupling term in the free energy rather than a competition 
between polar and rotational instabilities. 

We summarize the results up to this point by noting 
that the octahedral rotations and tilts in rhombohedral 
and c-oriented orthorhombic perovskites behave somewhat 
similarly in that compressive strain enhances octahedral 
rotations about the [001]-direction and tensile strain fa- 
vors a [110]-rotation axis with corresponding increases 
in the tilt angles. An exception occurs in etfr-oriented 
orthorhombic perovskites: the magnitude of both angles 
decreases with increasing tensile strain. In addition, we 
have seen in CaTiOs the first example of strain inducing 
a ferroelectric ground state. In the next section, we re- 
view further examples of such strong strain behavior in 
materials containing ferroelectrically active ions such as 
Ti 4 +. 
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C. Oxides with ferroelectrically-active ions and 
magnetic propensities 

Motivated by the observation in the previous section 
that strain can induce ferroelectricity in CaTi03, we next 
review additional examples that illustrate the interaction 
between strain and ferroelectricity. (For a thorough re- 
view see Ref. 11180 First we describe the behavior of 
diamagnetic SrTiC>3 and ferromagnetic EuTiC>3 - two 
additional cases which are not ferroelectric in the bulk, 
but in which the ferroelectrically active Ti 4+ ions are 
"activated" by strain.} ^ 9 * 120 * Next, we show how homoepi- 
taxial first-principles studies discovered that perovskites 
containing nominally Jahn- Teller inactive Mn 3+ could be 
coaxed to undergo ferroelectric displacements through epi- 
taxial strain constraints. And finally, we describe the case 
of BiFeC>3, which is already ferroelectric and magnetic in 
its bulk ground state, but where strain induces an u nusu al 
phase coexistence between two structural variantsP^ 

SrTiC>3 is not ferroelectric, but is an excellent 
dielectricP^I i n which the transition to a ferroelectric s tate 
is believed to be suppressed by quantum fluctuations^^ 
Instead, the ground state is tetragonal with an a°a°c _ tilt 
pattern of the oxygen octahedra. First-principles calcula- 
tions of bulk SrTi03 have shown that these antiferrodis- 
tortive rotations compete with an d ha ve a tendency to 
suppress the ferro electric inst ability P^J Early phenomeno- 
logical studies^ 2 ^! suggested that ferroelectric polarization 
could be obtained, and its orientation controlled, by ap- 
propriate strain conditions. 

Subsequent first-principles calculations^ 1 ^ were consis- 
tent with the phenomenological results. As in the rhom- 
bohedral perovskites, the DFT-based calculations find a 
change in the octahedral rotation axis from [001]— ► [110] 
on going from a compressive to tensile strain state. In ad- 
dition, polar displacements activated by epitaxial strain 
are found, with a [110] orientation of polarization fa- 
vored by tensile strain, as in the case of CaTiCv In 
contrast to the CaTiOs case, however, polarization is 
also induced by compressive strain, this time along the 
[001]-direction. Only between -0.4% and +0.2% strain 
is a paraelectric ground state found; in this region the 
out-of-phase rotations of oxygen octahedra dominated 
the structure. For larger compressive (or tensile) strains, 
a ferroelectric polarization is induced in the presence of 
these octahedral rotations. While ferroelectric hysteresis 
loops in strained SrTiOs have not yet been measured di- 
rectly, divergence of the dielectric constant, indicative of 
a ferroelectric phase transition, has indeed been reported 
at room temperature.^ 

The ability of first-principles calculations to identify 
critical strain regions where ferroelectric behavior is in- 
duced in common dielectrics has spawned numer ous ho- 
moepitaxial strain studies of magnetic dielectrics! 126 * 127 ! 
These ab initio searches for ferroelectricity in magnetic 
materials are motivated by the desire to identify new 
classes of magnetoelectric multiferroics, that is materials 
with simultaneous and coupled magnetic and ferroelectric 
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FIG. 11. Evolution of the lowest frequency infrared-active 
phonon with bi-axial strain for cubic EuTiOa with different 
magnetic configurations: spins aligned parallel (ferromagnetic) 
and antiparallel (G-type antiferromagnetic) . Strain states 
for which oj 2 < indicate a ferroelectric instability. Data 
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properties;- 

An example of a material studied with this motiva- 
tion is perovskite-structure EuTi03, which is isovalent 
with SrTi03 , but has the additional feature of magnetic 
/-electrons on the Eu 2+ ions. Bulk EuTiC>3 is reported 
to have the ideal cubic perovskite structure, with no 
ferroelectric polarization, and antiferromagnetic order- 
ing of the local Eu 2+ magnetic moments. Because the 
atomic and electronic structures of EuTiC>3 closely re- 
semble those of SrTi03, a similar strain-induced ferro- 
electricity should be anticipated in this case . Indeed, 
homoepitaxial first-principles calculations^ 2 ^ by Fennie 
and Rabe showed that application of ~1 % compressive 
strain is sufficient to cause the Ti 4+ cation to off-center in 
the direction perpendicular to the epitaxial plane — similar 
to the strain-induced ferroelectricity found in SrTi03 . As 
in SrTiC>3, the strain-induced ferroelectricity in EuTiC>3 
is understood to originate from the strong coupling be- 
tween strain-induced lattice deformations and the lowest 
frequency transverse optical modeP^l 

The authors of Ref. 11201 also found that under compres- 
sive strain the polar mode for a ferromagnetically ordered 
EuTiC>3 crystal is of lower energy (softer) than that of the 
antiferromagnetic spin arrangement at the same strain 
state (Fig. 
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This led Fennie and Rabe to suggest 
that strain could simultaneously modify both magnetic 
and electric ferroic orders. Recent calculations^ showed 
that an intriguing ferroelectric and ferromagnetic EuTiC>3 
phase should also be accessible under tensile strains larger 
than 0.6%. Using those ab initio guidelines, a stable mul- 
tiferroic phase, with mutual ferroic coexistence, has been 
subsequently realized experimentally.^ 

The perovskite-structure rare-earth and alkaline-earth 
manganites are of tremendous interest because of their 
rich structural, magnetic and electronic phase diagrams 
and magnetoresistive behavior. Incorporating ferroelec- 
tricity in the insulating members of the series would add 
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another desirable functionality. Under usual conditions, 
however, Mn 4+ and Mn 3+ cations do not undergo ferro- 
electric off-centering because the non-zero d-orbital occu- 
pation intr odu ces a large electronic penalty for off-centric 
distortionsP22lln this capacity, first-principles calculations 
have been used to explore circumstances under which such 
an off-centering might be induced. 

With increasing cation size, the bulk structures evolve 
from orthorhombic (Ca), with large octahedral rotations, 
to cubic (Sr), and finally a hexagonal structure (Ba), 
characterized by both corner- and edge-shared octahedra. 
CaMn03, although centrosymmetric, in fact has a ferro- 
electric instability in the cubic phase that is quenched by 
the a~a~c + octahedral rotations observed in the ground 
state structure Motivated by the sensitivity of the fer- 
roelectric mode to strain, Bhattacharjee and co-workers 
performed first-principles homoepitaxial strain calcula- 
tions for CaMnOa and identified that the competition 
between the rotational and ferroelectric instabilities fa- 
vors t he p olar structure for tensile strains greater than 
~ 2%1A2£I Here, a ferroelectric polarization develops in 
the ep itax ial plane, driven by displacements of the Mn 
cationJ^l that coexists with the competing a~a~c + octa- 
hedral rotation pattern. They also showed that compres- 
sive strain (up to 4%) does not stabilize the ferroelectric 
instability. Consistent with our earlier discussion, the 
authors of Ref. 11351 found that the frequency of the oc- 
tahedral rotational instabilities are less sensitive to the 
bi-axial strain than the frequency of the lowest polar 
phonon. 

Similar homoepitaxial strain calculations h ave been 
performed on perovskite-structured SrMn03p23 anc j 
BaMnOaP^In SrMnC>3, the Mn 4+ cations undergo ferro- 
electric off-centering under strain, even in the presence of 
octahedral rotations. Unlike CaMn03, however, SrMnC-3 
shows both out-of-plane (for compressive strains larger 
than 1.4%) and in-plane (tensile strains larger than 1%) 
polarizations. These critical strains are smaller than that 
for CaMnOs due to the larger Sr-cation which effectively 
produces an "internal" chemical strain on the lattice. In- 
terestingly, Lee and Rabe also report a large spin-phonon 
coupling, similar to that of EuTiC>3 for SrMn03 : strains 
greater than approximately ±3% induce a transition from 
an antiferromagnetic to ferromagnetic spin configuration 
on the Mn cations which are simultaneously displaced 
from the center of their oxygen coordinating octahedral 1 -^ 

In hypothetical perovskite-structure BaMn03, on the 
other hand, the larger Ba cation stabilizes the ferroelec- 
tric state even at its equilibrium volume — no strain is 
required!^ This occurs because the perovskite phase, 
which is metastable with respect to the denser hexago- 
nal structure, has such a large cell volume that the Mn 
cation becomes severely underbonded. Therefore, the Mn 
cation off-centers towards the edge of an octahedron in 
order to make two strong Mn-0 bonds. The authors of 
Ref. 08| also show that the ferroelectric perovskite struc- 
ture becomes the lowest energy phase at very large tensile 
strains, since under those elastic conditions, the denser 
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FIG. 12. (a) Evolution in the total energy for BiFeOs as 
a function of in-plane strain. The insets indicate the two 
structural variants — both with monoclinic symmetry — with 
the long (left) and short (right) c-axes that are accessible 
under large compressive and modest strain, respectively. In 
(b) the c/a axial ratio shows an abrupt discontinuity at the 
isosymmetric transition near -4.5% strain. The lattice strains 
corresponding to a number of commonly used oxide substrates 
are shown as dashed lines (see Table [IT] for a key to the labels). 
Data reproduced with permission from Ref. 11211 Copyright 
2009, American Association for the Advancement of Science. 



hexagonal phase is energetically unstable. Ferroelectric 
behavior has not yet been observed experimentally in 
any of these Mn-based compounds. We note that simi- 
lar volume-dependent ferroelectric instabilities have also 
been reported in chromate-based perovskitesP^ Substan- 
tial efforts are underway, however, to explore whether 
combinations of alkaline earth cations in manganite su- 
perlattices subjected to various bi-axial strain conditions 
can stabilize ferroelectric behavior. 

Finally for this section, we discuss the strain- 
dependence of the behavior in BiFeOs , in which the bulk 
ground state is already magnetic and ferroelectric, but 
which shows a strong evolution of the ferroelectric be- 
havior with strain. Bulk BiFeOs has the rhombohedral 
i?3c structure, which consists of antiferrodistortive octahe- 
dral rotations (a - a~a - )around the [111] axis, similar to 
those of LaNi03 and LaAlC>3 , and an additional relative 
off-centering of anions and cations along the [111] direc- 
tion leading to a ferroelectric polarization along that axis. 
In Figure 12 we show the calculated total energy (up- 



per panel) and c/a ratio (lower panel) for homoepitaxial 
BiFeOs films as a function of in-plane biaxial strain (from 



Ref. 11211 While the symmetry of the system remains mon- 
oclinic Cc throughout, there is an isosymmetric phase 
transition at ~4% compressive strain which is character- 
ized by an abrupt change in c/a ratio, and a change in the 
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coordination environment of the Fe from [6]-coordinated 
octahedral to [5]-coordinated square pyramidal P^ 21 The 
transition is accompanied by a re-orientation and enhance- 
ment of the ferroelectric polarization from ~ 90 /iC/cm 2 
about an axis close to [111] to 150 /xC/cm 2 about an axis 
close to [001]. There is also a change in the octahedral tilt 
pattern from a pattern that is strongly reminiscent of the 
bulk-like a~a~a~ tilt pattern to rotations mainl y ab out 
an axis that is close to the out-of-plane directionP^I We 
note, however, that there are multiple structurally unique 
states fRef. I140|) . which are less than 100 meV/f.u. higher 
in energy that compete with this isosymmetric phase tran- 
sition. If any of these phases with different symmetry 
occurs experimentally, the transition is of course no longer 
isosymmetric (see Refs. [TITI and IT4"IZ|) . While the strain 
value at which this transition is predicted to occur might 
be expected to be too large for experimental realization, 
in fact, thin films of BiFeC>3 grown on YAIO3 (with a 
5% lattice mismatch) do form in the large c/a structure. 
Even more intriguingly, films on LaAlOs, which h as a 
lattice constant corresponding to the cross-over point 
show a coexistence of the two phases with a so-called 
self-morphotropic phase boundary b etwe en them that can 
be manipulated by an electric field 

In summary, we have seen in these examples, first that 
imposing coherence with a substrate removes some of the 
symmetry (diad, triad, or tetrad) axes (Table |l| about 
which the octahedra rotate. Subsequently, biaxial strain 
modifies the rotation patterns by altering the magnitude 
of the rotation angles about these axes, and in the ex- 
treme case de-activating or activating new tilt patterns. 
In addition, bi-axial strain deforms the BOq octahedra 
by elongation or compression of the B-0 bond lengths. 
The structural distortion that dominates depends on the 
compressibility of the B-0 bonds and tendency for the 
octahedra to rotate as gleaned from temperatu re and 
pressure experiments on a range of perovskites! 144 * 14 ^ 
The changes in symmetry, bond angles and bond lengths 
in turn can have profound effects on the properties of 
the films. For example, responses that are prohibited 
by symmetry in the bulk may become allowed, changes 
in bandwidths can lead to drastically different electronic 
and optical properties, and changes in exchange inter- 
actions can change magnetic properties. The LaAlC>3 
and LaNiC>3 examples illustrated that, even in ostensible 
simple materials, the strain response is in general complex 
and simple models such as treating the oxygen octahedra 
as rigid units, or in the opposite extreme ignoring the 
response of the octahedral rotations may not be reliable. 
Materials that are proximal to ferroelectric instabilities, or 
that are already ferroelectric in their bulk ground states, 
show even more complex strain responses. Indeed, simple 
guidelines for the change in structure with strain are not 
yet available, and until a larger database is established 
we recommend full first-principles calculations with relax- 
ation of all variables rather than models or intuition for 
predicting the structural response of thin film oxides to 
strain. 



V. RESULTS: HETEROEPITAXIAL STRAIN. 

It is often asserted in the literature that the symmetry 
and structure of a substrate imprint across an interface 
so that a coherently grown film is affected not only by 
the substrate lattice constant but also by the details of 
its structure. In this section, we review the results of 
electronic structure calculations that have been designed 
to test this hypothesis, with a particular focus on the 
propagation of tilt patterns of oxygen octahedra across 
interfaces. We emphasize again that, in an electronic 
structure calculation, the effects of the presence of an 
interface can be studied independently from the effects of 
strain by comparing the outcomes of homo- and hetero- 
epitaxially strained systems. This provides valuable in- 
formation, which is difficult, if not impossible to obtain 
experiment ally. 

a. SrFe0 3 /SrTi0 3 . We begin with a study of 
a "model" heterostructure system: SrFe0 3 /SrTi0 3 
(SFO/STO), which is a good prototype system to evaluate 
the interplay of octahedral rotations across a heteroin- 
terface for a number of reasons: First, both materials 
are formed from neutral (001) planes of ^40 or BO2 
ions, therefore, the heteros truc ture avoids the complica- 
tion of a polar discontinuity^ 4 ^ (and in turn polar distor- 
tions) at the heterointerface between SrFeC>3 and SrTiC>3. 
Bulk SFO is metallic with p-type conductivity} 1 ^ and 
is proximal to multiple instabilities: it manifests a long- 
wavelength spin-density wave, but neither Jahn- Teller 
distorts nor charge orders, even though both possibili- 
ties are suggested by the high-spin d A chemistry of the 
Fe 4+ ion. SrTiC>3 is a highly polarizable dielectric, whi ch 

T48l 



can couple to electronic or structural distortions 1 ^ 1 in 
the SFO layer. The band alignment across the heteroin- 
terface is also Schottky-like within the LSDA, and the 
interface does not suffer from pathologies associated with 
the DFT underestimation of the band gapP^ Finally, for 
the study of octahedral rotations, this heterostructure 
is ideal because the bulk compounds exhibit simple oxy- 
gen octahedral tilt patterns: SrFeC>3 has the ideal cubic 
Pm3m perovskite structure (a°a°a° ti lt p attern) down to 
the lowest temperature studied (~4 KjE>2] anc j th e ground 
state 1 4/ mem phase of SrTiOs - which is a widely used 
substrate - has a single octahedral instability with re- 
spect to the cubic phaseP^ Below ^105 K it exhibits the 
a°a°c~ tilt pattern. 

In Ref. IBOl the authors investigated the effect 
of heterostructure periodicity in both symmetric 
(SrTi03)„/(SrFe03)„, n = 1...5, and asymmetric 
(SrTi0 3 )n/(SrFe0 3 ) m , n = 1...3,m = 1...3 superlat- 
tices. It was found that the octahedral rotations from 
the SrTiC>3 substrate propagate into the first two interfa- 
cial SrFeC>3 layers (Fig. 13), regardless of the number of 
SrTiC>3 layers. The rotational tendencies of the SrTiC>3 
layers are imprinted into the SrFeC>3 even in the (1,1) 
heterostructure. 

In heterostructures with ultra-thin (one- or two-layer 
thick) SrFeC>3 layers, these substrate-induced tiltings 
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FIG. 13. The layer-by-layer resolved octahedral rotation angles (8) for the(n, m) asymmetric (a) and (n, in = n) symmetric (b) 
(SrTi03)n/(SrFeOs) m heterostructures. The magnitude of the a°a°c- tilt pattern rapidly decreases away from the interfacial 
layers as shown schematically in (c) for the (SrFeOs)3/(SrTi03)3 heterostructure. Here, the magnitudes of the octahedral 
rotations about the axis (light, red line) perpendicular to the interface (bold line) are indicated by the length of the arrows. 
After Ref. 1601 Copyright 2010, American Physical Society. 



combined with the quantum confinement induce addi- 
tional instabilities - charge-ordering and/or Jahn- Teller 
distortions - that are not observed in bulk SrFe03. The 
authors® point out that the octahedral rotations are 
different for each electronic instability: charge ordering 
prefers the a°a°c~ tilt pattern while the Jahn- Teller dis- 
tortions is found to coexist with the a~a~c + tilt pattern. 
This ab initio result is consistent with a recent group theo- 
retical analysis of octahedral rotation s and electronically- 
driven structural distortions! 152 * 15 ^ It is worth noting, 
however, that while specific rotational patterns occur 
with each electronic instability it remains unclear whether 
the rotational pattern induces the electronic instability 
or vice versa. This merits additional study since it offers 
a possible route to controlling electronic phases through 
octahedral rot atio ns. (A similar suggestion has also been 
made recentljff^ for LaTiC>3 monolayers embedde d in 
SrTiC>3 and at manganite/titanate heterointerfacesl^J) 
In these highly confined ferrate heterostructures, the 
Jahn- Teller and charge orderings are accompanied by 
metal-insulator transitions in the nominally bulk metal- 
lic SrFeC>3 layer liSSl Corresponding homoepitaxial strain 
calculations show that the octahedral and electronic lat- 
tice instabilities are not induced in SrFe03 using bi-axial 
strain alone, indicating that substrate coherency and con- 
finement play a critical role in determining the structure 
and properties in these heterostructures. 

b. Manganite Superlattices. A more complicated 
model system is provided by LaMnOs/SrMnOs super- 
lattices, which combine magnetism with orthorhombic 
(a~a~c + ) and cubic (a°a°c°) symmetries (tilt patterns). 
While both constituents are antiferromagnetic insulators, 
there is additional electronic complexity introduced by 
the d 4 Mn 3+ ion in LaMn03 which has a tendency to 
Jahn- Teller distortion, and by the polar discontinui ty — 
LaMn03 has char ged (00 1) layers - at t he in terfaceP^ 
Both experimentaP^HSni anc \ theoretical^^ studies of 



these superlattices have focused on how the epitaxial 
strain coupling between the spin and orbital degrees of 
freedom at the heterointerfaces influences the macroscopic 
properties. The magnetism and orbital ordering are ex- 
pected to be highly sensitive to the strain condition at the 
interface, since changes in the bond angles and lengths 
will alter the preferred exchange mechanism and thus the 
flavor of orbital ordering. In addition, possible metallic- 
ity due to the polar discontinuity or interfacial mixing 
could change the dominant interactions from super- to 
double-exchange. 

Experimentally it has been found that the orbital degree 
of fr eedo m is indeed strongly modulated by the strain 
stateP^ In compression a C-type insulating antiferro- 
magnetic state is stabilized, while in contrast, tensile 
strain produces an A-type conducting interface. For the 
lattice-matched case, ferromagnetic order is observed, con- 
sistent with a disordered orbital state. First principles- 
calculations of layer-resolved band structures'^] indicate 
a crystal field degeneracy splitting induced by the inter- 
face strain that supports this interpretation of the above 
experiments; in addition a spin-polarized electron gas 
is calculated to occur at the interface in larger per iod 
manganite superlattices due to polar mismatch effects!^^ 
However, we emphasize that the calculations of Refs. 11621 
and 11631 did not allow for the presence of rotations or 
tilts of the oxygen octahedra, which we have seen can 
drastically alter the physics. In fact, modulations in the 
octahedral rotations in such manganite superlattices have 
been shown exp erimentally to alter the magnetic ordering 
temperature^^ — effects of local octahedral distortions in 
artificial structured materials are thus an obvious area 
for future study. 

There have been a number of recent experimental mea- 
surements of octahedral rotations across interfaces, that 
are broadly consistent with the picture that is emerg- 
ing from the first-principles calculations. For example 
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FIG. 14. In ultra-short period (SrTi0a)i/(PbTi03)i superlattices an unusual combination of out-of-phase (a) and in-phase (b) 
rotations combine about the axis perpendicular to the interface to support a ferroelectric polarization along the same direction 
(c). Here, the arrows indicate the rotation directions of the octahedra in each layer and their magnitude corresponds to their 
relative amplitude. After Ref. 3j Copyright 2008, Nature Publishing Group. 



real space mapping^ of the octahedral rotations across 
the La . 7 Sr .3MnO 3 /BiFeO3 (LSMO/BFO) interface with 
scanning transmission electron microscopy show the rota- 
tions propagating across the interface although modulated 
in magnitude from their bulk values to avoid energetically 
costly frustrations of the rotations across the heteroint- 
erface. And in (LaNi03)„/(SrMn03) m superlattices, the 
penetration length of the rotations across the interface has 
been shown to depend on the distance between different 
perovskite layers composed of large and small rotation 
anglesPH 

c. More exotic behaviors - some consequences of sym- 
metry lowering. Finally, we briefly mention some exotic 
behaviors in which the symmetry lowering associated with 
the presence of a heterointerface allows new properties 
to develop that are genuine properties of the interfacial 
system rather than either bulk parent compound. 

First-principles calculations were recently used to 
demonstrate a novel kind of improper ferroelectricity - in 
which the primary order parameter for the phase transi- 
tion is not the ferroelectric polarization - 1/1 period super- 
lattices of ferroelectric/paraelectric PbTiOs/SrTiOaP In 
this case, tetragonal ferroelectric PbTi03 (which does not 
have octahedral rotations) is combined with antiferrodis- 
tortive SrTi03 and an enhanced polarization is obtained; 
the conventional bulk description of ferroelectricity would 
suggest that interruption of the cooperative Ti displac- 
ments in the PbTi03 layers by paraelectric SrTi03 should 
attenuate the macroscopic polarization. The authors of 
Ref. |3| showed that because of the competition among the 
structural instabilities at the heterointerface, an unusual 
antiferrodistortive rotation of the oxygen octahedra about 
an axis perpendicular to the interface is stabilized. This 
is symmetry-compatible with a ferroelectric polarization 
along that direction (Fig. [M| . Bousquet and co-workers 
demonstrated theoretically that the enhanced polariza- 
tion is driven by the specific octahedral rotation pattern 
that is present at the heterointerface but absent in the 
bulk constituents. Motivated by the predictions, thin 
films of PbTiOs/SrTiOs were grown and a very large 



and temperature-independent dielectric constant, typical 
of improper ferroelectrics but unusual for conventional 
ferroelectrics, was measured P 

A similar enhancement in interfacial ferroelectric po- 
larization was also predicted using heteroepitaxial DFT- 
based calculations for asymmetric and symmetric superlat- 
tices of paraelectric CaTi03 and ferroelectric BaTiOsP^S 
Here, Wu et al. found that large TiOg octahedra rota- 
tions persist in all superlattices studied when the adjacent 
oxide layers are CaO, but the amplitude of the rotations is 
substantially reduced when the adjacent layers are BaO. 
This behavior is consistent with the bulk structure of 
CaTi0 3 (BaTi0 3 ) which has the a^a~c+ (a°a°c°) tilt 
pattern. With an increasing ratio of CaO to BaO layers, 
or vice versa, the rotation magnitudes approach their re- 
spective bulk values; however, in the ultra-short 1/1 limit 
the rotations are about half the size (4° - 6°) of those in 
bulk CaTiOaP^I Because the octahedral rotations and 
ferroelectric displacements compete with each other in 
bulk CaTi03, the suppressed octahedral rotations at the 
heterointerface between BaO and CaO layers allow for a 
larger polarization to develop in that layer, which in turn 
enhances the net polarization of the superlattice. Since 
these structural distortions strongly couple to bi-axial 
strain, they could be further enhanced by growth on a 
suitable substrate. 

In the previous two examples, the oxide layers at the 
heterointerface break inversion symmetry in the synthetic 
perovskite because the chemical and structure environ- 
ments in directions perpendicular to the interface are in- 
equivalent. Since the linear magnetoeletric effect can only 
be non-zero in the absence of time-reversal and space- 
inversion symmetries, heterointerfaces can be used to 
enable magnetoelectric response in otherwise centrosym- 
metric magnetic materials. In Ref. [168, first-principles 
calculations were used to demonstrate such a linear mag- 
netoelectric effect at the interfaces in SrRuOs/SrTiOs 
superlattices. The effect is symmetry prohibited in both 
parent compounds, but is allowed at the interface. The 
authors demonstrated that the magnetoelectric response 
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arises from a carrier-mediated mechanism, and should be 
a universal feature of the interface between a dielectric 
and a spin-polarized metal; it has subsequently b een con- 
firmed in a CoPd film immersed in an electrolyte^^ anc j 

all-solid- state fe rromagnetic (La,Sr)Mn03/Pb(Zr,Ti)C>3 
interface™!!] 

Similar calculations have now been performed on a 
range of ferroelectric/ferromagnetic metal interfaces and 
novel interfacial multiferroic behavior reported (see for ex- 
ample Refs. ll72lH75H t is important to note however, that 
calculations for ferroelectric/metal interfaces are fraught 
with technical difficulty, because the DFT underestima- 
tion of the electronic band gap in the insulator can lead 
to calculated ohmic contacts in situations where a Schot- 
tky barrier occurs experimentally. As a result, spurious 
real-space charge transfer occurs and this can obscure 
the intrinsic behavior of perovskite heterointerfaces with 
competing structural and electronic instabilities. For a 
detailed discussion of the unphysical behaviors caused 
by this pathology, including many examples from the 
existing literature, see Ref. 11491 Methods such as the 
recently developed formalism for performing density func- 
tional calculations for capaci tors w ith constrained values 
of the dielectric displacement 65 ! 17 -^ go some way towards 
alleviating this problem. 

VI. OUTLOOK FOR RATIONAL OXIDE 
HETEROSTRUCTURE DESIGN 

In this review we described how first-principles calcula- 
tions based on density functional theory can be used to iso- 
late atomic and electronic structure changes in perovskite 
oxide thin films and heterointerfaces. We described ef- 
forts using the homoepitaxial and heteroepitaxial strain 
approaches to decouple the intrinsic contributions that 
epitaxial strain, changes in symmetry, and interface chem- 
istry play in determining the properties of oxide het- 
erostructures. 

Several new ideas for oxide heterostructure and thin 
film design emerge from this review: 

o The macroscopic properties of oxide heterostruc- 
tures can often not be simply predicted from con- 
sideration of the electronic structure of the bulk 
materials alone; the interfacial physical, electronic 
and magnetic structures in artificial geometries can 
be genuinely different from those of the parent bulk 
materials due to changes in symmetry- and size- 
dependent properties. 

o Bi-axial strain does more than just change bond 
lengths in perovskite thin films; it can couple to 
and/or alter the internal degrees of freedom. In 
particular, octahedral rotation patterns are modu- 
lated by strain in a fashion that is not immediately 
intuitive and can give rise to new electronic states. 

o The heterostructure ground state is often influenced 
by latent instabilities present in the bulk phases; 



the substrate-induced heteroepitaxial constraints 
then act to enhance any unstable modes and to 
re-normalize the low energy electronic structure. 

o Translational and point symmetry changes at a het- 
erointerface can lift bulk electronic degeneracies 
to promote new and allowed order parameter cou- 
plings, through for example, strain-induced crystal 
field splittings, which in turn strongly affect the 
macroscopic behavior. 

Future research directions. Finally we outline some 
on-going and new research directions and pressing open 
questions in the field that we find of particular interest. 

1. Mechanisms for strain accommodation. We have 
seen that structurally similar perovskite oxides be- 
have differently under epitaxial strain both from 
each other, and from their bulk counterparts un- 
der hydrostatic pressure. While much of the elastic 
strain accommodation in heteroepitaxy is accom- 
modated through changes in both the bond lengths 
and octahedral rotation patterns, it is still unclear 
why some oxides show relatively larger bond length 
changes, while others undergo larger changes in the 
octahedral rotation and tilt patterns. The answer 
likely is found in the different compressibility of cer- 
tain transition metal - oxygen bonds; future work 
should attempt to quantify this by surveying a va- 
riety of insulating and metallic perovskite oxides, 
and analyzing the stiffnesses of octahedral rotations 
and bond length distortions. By collecting these 
data, it may be possible to build a set of design 
rules governing the tendency for certain classes of 
perovskite oxides to undergo different atomic dis- 
placement patterns with bi-axial strain. 

2. Octahedral texturing. We have seen many exam- 
ples in which bi-axial strain strongly modifies the 
octahedral rotation patterns in thin films by alter- 
ing the symmetry axes about which the octahedra 
rotate. In all cases, when the in-plane lattice pa- 
rameters are equal and the transition metal nearest 
neighbor distances are the same, then the rotations 
around the x- and y-axes are equal (or equivalently, 
the net rotation is about the [110]-direction). How- 
ever, when the two in-plane lattice parameters differ, 
this degeneracy is lifted and different directions of 
rotation are adopted.^ This geometric constraint 
suggests that the substrate miscut angle - which 
can change the effective in-plane a and b lattice 
parameters adopted by a coherent film - could be 
a useful parameter to control the octahedral rota- 
tions and the preferred rotational easy axis. Such 
crystallographic substrate-tailoring could make it 
possible to achieve monodomain samples, or obtain 
a specific number of antiphase domain boundaries. 
It has recently been demonstrated experimentally 
that substrate orientation can be used to modify the 
orbital ordering patterns in manganite compounds: 
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a "striped" phase can be transformed to a "checker- 
board" phase by switching from growth on a (001) 
to (110) terminated surfaceP^ The two surface ter- 
minations in this case reflect the extreme situation 
of a substrate miscut angle, and we propose that 
by slowly tuning that angle, both orbital polariza- 
tion and octahedral rotation orientation might be 
controllable in related thin film oxides. 

3. Intrinsic defect profiles. An important ques- 
tion, which in principle is accessible through first- 
principles calculations, is whether the intrinsic de- 
fect profile of a film changes with strain or heteroepi- 
taxy. It is widely believed in the oxides community 
that the lattice constant of perovskite oxides de- 
pends on the concentration and type of intrinsic 
defects Indeed, accurate measurements of lattice 
constants have sometim es been used to infer oxygen 
vacancy concentrations P*1!I2I Therefore, it is likely 
that at large strain values, it may become more en- 
ergetically favorable to accommodate changes in lat- 
tice constants through changes in the defect profile 
rather than in changes of the bond lengths and tilt 
angles — a chemical strain rel axation mechanisirP^I 
consistent with Vegard's lawf^MHH j n addition to 
being of fundamental interest, this issue is of pro- 
found technological importance since the properties 
of oxide films often depend sensitively on defect 
concentrations. While in an experiment, multiple 
defects - for example oxygen vacancies and cation 
non-stoichiometry - usually occur simultaneously, 
in calculations it is possible to separately evaluate 
the effect of each individual defect on the lattice 
parameters of the system. A comprehensive series of 
such calculations would be helpful in identifying the 
intrinsic changes in lattice parameters with defect 
profile, and in turn likely changes in defect profile 
with strain. In practice, however, calculations for 
realistic defect profiles require large supercells, and 
accurate calculations of stresses induced by the in- 
troduction of defects require large energy cut-offs. 
Therefore such a study is a formidable task for fu- 
ture work. 



4. Technical issues. It is a well-known problem within 
density functional theory that standard exchange- 
correlation functionals such as the local density 
approximation have errors in their calculation of 
lattice constants of up to a few percent. These er- 
rors can in turn have rather drastic consequences 
on the properties. For example, the prototypical fer- 
roelectric BaTiC>3 is paraelectric at the theoretical 
local density approximation (LDA) lattice constant, 
while paraelectric SrTiC>3 is f erro electric at the the- 
oretical GGA lattice constant!-^ In standard "bulk" 
calculations, this is often circumvented by working 
at the experimentally measured lattice constants. 
This is not a possibility, however, when strain is to 
be used as a variable in a calculation, as the experi- 
mental lattice parameters have a theoretical strain 
associated with them, and the only well-defined 
zero-strain reference state is the calculated struc- 
ture. Development of exchange-correlation func- 
tionals that yield accurate bulk lattice parameters 
and structures, and testing of their behavior in 
strained systems, is therefore a crucial direction for 
future research. Here, the recently introd uced GGA 
exchange-correlation functional PBEsol, 18 which 
is biased towards more accurately reproducing sur- 
face energies and lattice constants than its PBE 
(Perdew-Burke-Ernzerhof ) pr edecessor seems to 
be particularly promisingPSSl 

In summary, first-principles studies have shown that 
prior conventional wisdom guiding TMO heterostructure 
design as simple two-component composites should be re- 
evaluated. Instead, understanding the electronic phases at 
oxide heterointerfaces requires self-consistent treatment of 
the electronic and atomic degrees of freedom of both con- 
stituents on an equal footing. We have seen many exam- 
ples in which, because of the many competing electronic 
and structural degrees of freedom in perovskite oxides, 
the physical properties found in oxide-oxide heterostruc- 
tures are highly susceptible to subtle changes in elastic 
strain and dimensionality. This critical understanding of 
how correlated electron and emergent behavior develops 
from changes in local structure and artificial geometries 
is essential to engineering their functionality. 
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